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Abstract 

We show that the presence of residual Bose-Einstein correlations 
may affect the non-resonant contribution of hadronic charm decays 
where two identical pions appear in the final state. The distortion of 
the phase-space of the reaction would be visible in the Dalitz plot. 
The decay Z)+ K^tt^tt^ is discussed but results can be generalized 
to any decay with identical bosons. 
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1 Introduction 



Hadronic decays of D mesons are an important source of information. The 
study of pseudoscalar-pseudoscalar and pseudoscalar- vector decays of D mesons, 
for example, shed hght on heavy quark decay and hadronization. 
The hfetime difference between the and the D'^ mesons as well as the 
role of final state interactions may be understood by carefully studying their 
hadronic decays. 

The fractions of resonant (pseudoscalar-pseudoscalar or pseudoscalar- 
vector) and non resonant three-body contributions to each decay mode are 
determined by fitting the Dalitz plot distributions to a coherent sum of am- 
plitudes. 

The three-body amplitude is normally assumed to be constant over the Dalitz 
plot but, decays like —>■ K'h^tt^ contain identical charged pions in the 
final state and their interference produce a non-constant behaviour of the 
phase-space p. 

The first evidence for a non-uniform population on the Dalitz plot for 
K'-K^Ti^ was shown in ref. ^j. 

In an attempt is made to account for the interference by adding an ad 
hoc term to the likelihood function. The best fit is obtained with a Bose 
symmetric term of the form B^body =\ ^i-jri ~ ''^I:7r2 I- However, they found 
that this form does not seem to adequately describe the physics. 
In a recent analysis j3| the amplitudes of the fit were again unsuccessfully 
Bose-symmetrized in an attempt to describe the non-resonant contribution. 

Recently Bediaga et al. [T studied the problem and found that the non- 
resonant contribution should actually not be constant over the phase-space. 
They used a method based on factorization and an effective hamiltonian for 
the partonic interaction to estimate the distortion of the Dalitz plot. 

Here we will try to address the problem from the perspective of Bose- 
Einstein interference among the identical bosons in the final state and the 
residual correlation produced. This may help to understand the Dalitz plot 
of K'-K^-K^ in particular, and the effect in any decay with several 

bosons in general. 

2 Bose-Einstein correlations 

The study of Bose-Einstein correlations (BEC) in connection with other phe- 
nomena [5| in high energy physics reactions became important recently when 
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it was realized that such correlations may affect measurements of the stan- 
dard model parameters [HIIZI- 

Charm mesons travel long enough to neglect the effect of interference 
among their decay products and the pions produced in the reaction. Effects 
of the kind as described in [Sj will not be considered here. We rather will 
regard the decay itself as a particle production process in which the interfer- 
ence may arise. 

Fig. 1(a) shows the decay with a bubble representing the region in space- 
time where particles are produced. Decays where particles are produced via 
fragmentation are not completely coherent. The incoherence present in the 
decay gives rise to the Bose -Einstein interference. 

The Bose-Einstein correlations are commonly described in terms of a two 
particle correlation function: 



where P{pi,p2) is the joint probability amplitude for the emission of two 
bosons with momenta pi and p2, and P{pi) and P{p2) are the single produc- 
tion probabilities. 

The Bose-Einstein correlation among identical mesons has been used to 
probe the space-time structure of the intermediate state right before hadrons 
appear [HI IH! in high energy and nuclear collisions. 

Real data are analysed in terms of the ratio of distributions of like charged 
pairs (identical bosons) to that of an uncorrelated sample 



The uncorrelated sample Co(pi,p2) with all the features of C{pi,p2) except 
for the Bose-Einstein effect may be obtained from the distribution of unlike 
charged pions or by using an event mixing technique. 

One parametrizes the effect assuming a set of point-like sources emmiting 
bosons. These point like sources are distributed according to a density p(r). 
The correlation function can then be written as, 



Rbe{Pi,P2) 



(1) 







(3) 
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where pi,p2 are the momenta of the two bosons, i/jbe represents the Bose- 
Einstein symmetrized wave function of the bosons system and Jy p{r)d^r = 
1. Taking plane waves to describe the bosons one obtains the correlation 
function for an incoherent source: 

RBE{p[,p'2) = l+\Hp{r))\\ (4) 
where J-'{p) represents the Fourier transform of the density function p{r). 

In order to describe the quantum interference during the fragmentation in 
high energy reactions, phenomenological parametrizations of the effect have 
been proposed. For a recent review see pUj . 

One of the most commonly used parametrizations is given by: 

^BE(pi,P2) = l + Ae-«'^ (5) 

where is the Lorentz invariant, = — (pi — P2), which can be written 
also as, = rn? — 4m^, where m is the invariant mass of the two pions and 
uIt^ = O.lSOGey the mass of the pion. The parameter A lies between and 
1 and reflects the degree of coherence in the pion production. The radius of 
the pions source will be given by r = hc^/J3[fm\. 

The presence of Bose-Einstein correlations will modify not only the in- 
variant mass spectrum of like charged but also that of unlike charged pions. 
This reflection of BEG known as residual correlation has been studied since 
long time ago to make sure that the reference sample of unlike charged pions 
used to subtract the effect from the like charged pions spectra is clean. 
In particle production processes with high multiplicity, residual correlations 
easily dissolve. Altough some studies claim that this is not the case and 
that using unlike charged pions as a reference sample is questionable. 
In a particle production process where only three particles are yield {e.g. 
K^-K^Ti^) residual correlations are expected to be larger. 

3 Simulation of Bose-Einstein correlations 

In this letter we want to estimate the effects of BEG on the phase space of a 
three body decay. We take the particular case of — > K^ti^ti^ and look at 
it, as a particle production process with only two identical bosons in the final 
state. The reduced number of pions makes easy the incorporation of BEG 
in a Monte Garlo simulation. We will take the approach used in jS] where 
the BEG are simply simulated by weighting each event. The MG generator 
consists of a three body decay simulation with the appropiate masses of the 
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decay products. In j^] each event is weighted according to: 

= [](! + Ae-'^^'O, (6) 

where the product would be taken over all pairs of like charged pions. 
In our case the product is reduced to one pair of like charged pions 

iy = l + Ae-^«'. (7) 

We did simulations with different values for A and f3. Fig. 2 shows the 
squared invariant mass distributions both for K^n^ and 7r"'"7r"'" before and 
after simulation of BEC. The (5^(vr+7r+) distribution can be normalized using 
its distribution in the abscence of BEC. Fig. 3 shows the correlation function 
Rbe{t''~^'^~^) obtained according to eq. (2) and the curve of eq. (5) with A = 1 
and P = AGeV'"^. One can see that one obtains from the simulation exactly 
what one puts in. 



4 Effects on the Dalitz plot 

The effect of BEC on the Dalitz plot in some cases would be a straightfor- 
ward manifestation of the interference among identical pions. An example of 
this is the Dalitz plot of K'tt^h^ where the invariant mass squared 

of the two pions is plotted agains the invariant mass squared of the kaon and 
one of the pions, i.e. m^lTrn) vs. rn?{KTx). In this case the BEC will be 
visible directly in m?{mi). The invariant mass spectra of m^{K'K), however, 
will show residual effects. 

The Dalitz plot of rn?{KTs:i) vs. rn?{KTi2) will show only residual effects. 

Fig. 4 shows the Dalitz plots for the K^n^n^ decay once the 

BEC has been incorporated. The values A = 1 and P = AGeV~'^ were used 
to produce the distributions shown. Increasing (3 from 4:GeV~^ to lOOGeV^"^ 
would make the long bands in Fig. 4 narrower. As mentioned above, the 
physical meaning of /? = 4, 9, 16, 25, lOOGeV"^ is given in terms of the pions 
source radius r ^ 0.4, 0.6, 0.8, 1.0, 2. 0/m respectively. We show the results 
using r = 0.4/m just because the effect is more visible. Of course the real 
value would be extracted from fits to experimental data. 
A better fit of the non-resonant part of the decay would be possible once r 
has been extracted from data. 
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As in the case of the parameter /3, we used A = 1 because this pro- 
duces the biggest possible effect. The value A = 1 indicates a totally chaotic 
production in the decay. Of course a particle decay is not necessarily a com- 
pletely chaotic process. In fact, if fragmentation would not take place during 
the decay, one may expect a completely coherent process in which A = 
and Bose-Einstein interference would not be present. Fragmentation in the 
decay introduces some degree of incoherence giving to A a value between 
and 1. The exact value would be obtained fitting the correlation function 
as in the case of the source radius. On the other hand, there may be pro- 
duction mechanisms other than fragmentation that introduce some degree of 
incoherence. The study of BEG in particular decays may help to disentangle 
the production process. 

The hadronic decay —>■ K^-k^tt^ seems to be dominated by the 
non resonant contribution JJj and is therefore an interesting laboratory to 
study the role of BEG in hadronic decays in general. A resonant decay like 
— s> K*^Tx^ corresponds to a situation illustrated in Fig. 1(b). The in- 
terference between the pions in this case would, according to Bowler jT2], 
reduce and narrow the strength of the correlation. 

In hadronic decays where the resonant contribution is significant such effects 
must be taken into account. 

An interesting aspect of the fact that BEG correlations are observable in 
hadronic decays of charm and eventually beauty mesons is that there must 
be some relationship between the parametrization of p{r) for the pions source 
and the D meson form factors. Such relationship can be examined and the 
study of BEG may be a tool to learn more about the decay and fragmentation 
mechanism. 
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Figure Captions 

Fig. 1: Picture of (a) the non- resonant part of the decay in which the partonic 
processes are overlooked. The bubble described by p(r), parametrizes 
the space time of the production mechanism. In (b) a resonant decay 
is illustrated. 

Fig. 2: Squared invariant mass distributions for K^tt^ and tt^tt^ with (his- 
togram) and without (crosses) BEC. 

Fig. 3: Correlation function Rbe obtained normalizing the Q^^n^n^) distri- 
bution in the presence of BEC with the (^^(Tr+Tr"'") distribution without 
BEC. The curve corresponds to eq. (5) with A = 1 and P = AGeV~'^. 

Fig. 4: Dalitz plots for the — > K'^k^t:^ decay once the BEC has been 
incorporated. 
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Figure 4. 
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